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The theory presented here, cosmological general relativity, uses a Rie- 
mannian four-dimensional presentation of gravitation in which the coor- 
dinates are those of Hubble, i.e. distances and velocity rather than the 
traditional space and time. We solve the field equations and show that 
there are three possibilities for the Universe to expand. The theory de- 
scribes the Universe as having a three-phase evolution with a decelerating 
expansion, followed by a constant and an accelerating expansion, and it 
predicts that the Universe is now in the latter phase. It is shown, as- 
suming Qrn = 0.245, that the time at which the Universe goes over from 
a decelerating to an accelerating expansion, i.e., the constant-expansion 
phase, occurs at 8.5 Gyr ago. Also, at that time the cosmic radiation 
temperature was 146K. Recent observations of distant supernovae imply, 
in defiance of expectations, that the Universe's growth is accelerating, 
contrary to what has always been assumed, that the expansion is slow- 
ing down due to gravity. Our theory confirms these recent experimental 
results by showing that the Universe now is definitely in a stage of accel- 
erating expansion. The theory predicts also that now there is a positive 
pressure, p — O-OSig/cm^ , in the Universe. Although the theory has no 
cosmological constant, we extract from it its equivalence and show that 
A = 1.934 X lO-^'^s"^ This value of A is in excellent agreement with the 
measurements obtained by the High-Z Supernova Team and the Supernova 
Cosmology Project. It is also shown that the three-dimensional space of 
the Universe is Euclidean, as the Boomerang experiment shows. Com- 
parison with general relativity theory is finally made and it is shown that 
the classical experiments as well as the gravitational radiation prediction 
follow from the present theory, too. 



PACS numbers: 04.50.+li, 11.25.Mj., 11.27.+d, 98.80Cq 



*Email: car melim@bgumail.bgu.ac.il 



1 



1 Preliminaries 



As in classical general relativity we start our discussion in flat spacevelocity 

which will then be generalized to curved space. 

The flat-spacevelocity cosmological metric is given by 

ds^ = T^dv^ - {dx^ + dy^ + dz^) . (1) 

Here r is Hubble's time, the inverse of Hubble's constant, as given by mea- 
surements in the limit of zero distances and thus zero gravity. As such, r is a 
constant, in fact a universal constant (its numerical value is given in Section 8, 
T = 12.486Gyr). Its role in cosmology theory resembles that of c, the speed of 
light in vacuum, in ordinary special relativity. The velocity v is used here in the 
sense of cosmology, as in Hubble's law, and is usually not the time-derivative of 
the distance. 

The Universe expansion is obtained from the metric (1) as a null condition, 

ds = 0. Using spherical coordinates r, 0, (j) for the metric (1), and the fact that 
the Universe is spherically symmetric {dO = dtj) = 0), the null condition then 
yields dr/dv = r, or upon integration and using appropriate initial conditions, 
gives r = TV or v = H^r, i.e. the Hubble law in the zero-gravity limit. 

Based on the metric (1) a cosmological special relativity (CSR) was pre- 
sented in the text [1] (see Chapter 2). In this theory the receding velocities of 
galaxies and the distances between them in the Hubble expansion are united 
into a four-dimensional pseudo-Euclidean manifold, similarly to space and time 
in ordinary special relativity. The Hubble law is assumed and is written in 
an invariant way that enables one to derive a four-dimensional transformation 
which is similar to the Lorentz transformation. The parameter in the new trans- 
formation is the ratio between the cosmic time to r (in which the cosmic time 
is measured backward with respect to the present time). Accordingly, the new 
transformation relates physical quantities at different cosmic times in the limit 
of weak or negligible gravitation. 

The transformation between the four variables x, y, z, v and x', y', z', v' 
(assuming y' = y and z' = z) is given by 

,_x-t^ v-tx/r^ 

Equations (2) are the cosmological transformation and very much resemble the 
well-known Lorentz transformation. In CSR it is the relative cosmic time which 
takes the role of the relative velocity in Einstein's special relativity. The trans- 
formation (2) leaves invariant the Hubble time r, just as the Lorentz transfor- 
mation leaves invariant the speed of light in vacuum c. 



2 Cosmology in spacevelocity 

A cosmological general theory of relativity, suitable for the large-scale struc- 
ture of the Universe, was subsequently developed [2-5]. In the framework of 
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cosmological general relativity (CGR) gravitation is described by a curved four- 
dimensional Riemannian spacevelocity. CGR incorporates the Hubble constant 
T at the outset. The Hubble law is assumed in CGR as a fundamental law. 
CGR, in essence, extends Hubblc's law so as to incorporate gravitation in it; 
it is actually a distribution theory that relates distances and velocities between 
galaxies. The theory involves only measured quantities and it takes a picture of 
the Universe as it is at any moment. The following is a brief review of CGR as 
was originally given by the author in 1996 in Ref. 2. 

The foundations of any gravitational theory are based on the principle of 
equivalence and the principle of general covariance [6]. These two principles 
lead immediately to the realization that gravitation should be described by a 
four-dimensional curved spacetime, in our theory spacevelocity, and that the 
field equations and the equations of motion should be written in a generally 
covariant form. Hence these principles were adopted in CGR also. Use is made 
in a four- dimensional Riemannian manifold with a metric gr^i^ and a line element 
ds^ = g^j^dx^dx'^ . The difference from Einstein's general relativity is that our 
coordinates are: x^ is a velocitylike coordinate (rather than a timelike coordi- 
nate), thus = TV where r is the Hubble time in the zero-gravity limit and v 
the velocity. The coordinate x'^ = tv is the comparable to x'^ = ct where c is the 
speed of light and t is the time in ordinary general relativity. The other three 
coordinates x'^, k = 1,2,3, are spacelike, just as in general relativity theory. 

An immediate consequence of the above choice of coordinates is that the 
null condition ds = describes the expansion of the Universe in the curved 
spacevelocity (generalized Bubble's law with gravitation) as compared to the 
propagation of light in the curved spacetime in general relativity. This means 
one solves the field equations (to be given in the sequel) for the metric tensor, 
then from the null condition ds = one obtains immedialety the dependence of 
the relative distances between the galaxies on their relative velocities. 

As usual in gravitational theories, one equates geometry to physics. The first 
is expressed by means of a combination of the Ricci tensor and the Ricci scalar, 
and follows to be naturally either the Ricci tracc-frcc tensor or the Einstein 
tensor. The Ricci trace-free tensor does not fit gravitation in general, and 
the Einstein tensor is a natural candidate. The physical part is expressed by 
the energy-momentum tensor which now has a different physical meaning from 
that in Einstein's theory. More important, the coupling constant that relates 
geometry to physics is now also different. 

Accordingly the field equations are 

Gfj,iy = Rfj,iy —g^iyR - ^T^u, (3) 

exactly as in Einstein's theory, with k given by k = Sirk/r'^, (in general relativity 
it is given by SttG/c^), where k is given by fc = Gt'^ jf? , with G being Newton's 
gravitational constant, and t the Hubble constant time. When the equations of 
motion will be written in terms of velocity instead of time, the constant k will 
replace G. Using the above equations one then has k = SttG/c^t'^. 
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The energy-momentum tensor T'^'^ is constructed, along the hnes of general 
relativity theory, with the speed of light being replaced by the Hubble constant 
time. If p is the average mass density of the Universe, then it will be assumed 
that T^"^ = pu'^u" , where ?y/' = dx^^ /ds is the four- velocity. In general relativity 
theory one takes Tg = p. In Newtonian gravity one has the Poisson equation 
V^^ = AirGp. At points where p = one solves the vacuum Einstein field 
equations in general relativity and the Laplace equation V^(/) = in Newtonian 
gravity. In both theories a null (zero) solution is allowed as a trivial case. In 
cosmology, however, there exists no situation at which p can be zero because the 
Universe is filled with matter. In order to be able to have zero on the right-hand 
side of Eq. (3) one takes Tq not as equal to p, but to peff = p — Pc, where Pc 
is the critical mass density, a constant in CGR given by pc = 3/8itGt^, whose 
value is pc « 10~^®gf/cm^, a few hydrogen atoms per cubic meter. Accordingly 
one takes 

T^-- = PeffuV; Peff =p-pc (4) 

for the energy- momentum tensor. 

In the next sections we apply CGR to obtain the accelerating expanding 
Universe and related subjects. 



3 Gravitational field equations 

In the four-dimensional spacevelocity the spherically symmetric metric is given 

by 

ds'^ = T^dv^ - ei'dr^ - [de^ + sin^ ed(t>^) , (5) 

where p and R are functions of v and r alone, and comoving coordinates 
x'^ = {x^,x^,x^,x^) = {TV,r,9,(t)) have been used. With the above choice 
of coordinates, the zero-component of the geodesic equation becomes an iden- 
tity, and since r, 6 and (f) are constants along the geodesies, one has dx^ = ds 
and therefore 

= = (1,0,0,0). (6) 

The metric (5) shows that the area of the sphere r = constant is given by 
47ri?2 and that R should satisfy R' = dR/dr > 0. The possibility that R' = 
at a point rg is excluded since it would allow the lines r = constants at the 
neighboring points ro and ro + dr to coincide at ro, thus creating a caustic 
surface at which the comoving coordinates break down. 

As has been shown in the previous sections the Universe expands by the 
null condition ds = 0, and if the expansion is spherically symmetric one has 
dO = d(j) = 0. The metric (5) then yields 

T^rfv^ - e'^dr^ = 0, (7) 

thus 

? = re-^'''- (8) 
dv 
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This is the differential equation that determines the Universe expansion. In 
the following we solve the gravitational field equations in order to find out the 
function ^ (r.v). 

The gravitational field equations (3), written in the form 



(9) 



where 



PeffU^u^ + p (Uf^u^ - g^i,) . 



(10) 



with peff = p — and T = T^iyg^" , are now solved. Using Eq. (6) one finds 
that the only nonvanishing components of T^j/ are Too = r'^Peff, Tn = c ^rpe^", 
T22 = cr'^TpR^ and Tag = cr^rpR^ sin^ 6, and that T = r'^Peff - 3c" Vp. 

The only nonvanishing components of the Ricci tensor yield (dots and primes 
denote differentiation with respect to v and r, respectively), using Eq. (9), the 
following field equations: 



i?, 



00 



1 2 •• 

-Li- —R 
2^ R 



1,9 ^ / 2t o 1 \ 

= 2 V Peff + 3c TP) , 



i?oi = ^-^'A ~ -^R' = Oi 



1 



1 



R' 



1 



R 



R22 = RR+ ^RRp + + 1 - e"'^ (rr" - ^RR'p' + R'^^ 
R33 = sin^ 0R22 = sin^ e {r^Peff - c'^rp) . 



(Ua) 
(116) 

(11c) 

(lid) 
(lie) 



The field equations obtained for the components 00, 01, 11, and 22 (the 33 
component contributes no new information) are given by 



-/* - ^-R - = K [r^Peff + 3c Vp) 



27?' - R'fi = 0, 



1 



2 • 



/i + -/i^ + -Rp + e-^ —R'p' - —R 



R 



R 



R 



K {t^P, 



'eff 



-'rp) 



(12) 

(13) 
(14) 



1^- 1 



R'p' - 2 



5 



'eff 



(15) 



It is convenient to eliminate the term with the second velocity-derivative of /i 
from the above equations. This can easily be done, and combinations of Eqs. 
(12)-(15) then give the following set of three independent field equations: 



e'' {2RR + + l) - = - 
2R' - R'fi = 0, 



-^e'^R^p, 



+ 



i?2 



= peff, 

other equations being trivial combinations of (16)-(18). 



(16) 
(17) 

(18) 



4 Solution of the field equations 

The solution of Eq. (17) satisfying the condition ii' > is given by 

l + /(r)' 



(19) 



where / (r) is an arbitrary function of the coordinate r and satisfies the condition 
/ (r) + 1 > 0. Substituting (19) in the other two field equations (16) and (18) 
then gives 

2RR + R^ - f = -kc-WR^p, (20) 
^(2EE'-/') + i,(i?^-/)=.rVe//, (21) 

respectively. 

The simplest solution of the above two equations, which satisfies the condi- 
tion i?' = 1 > 0, is given by 

R = r. (22) 
Using Eq. (22) in Eqs. (20) and (21) gives 

f{r) = KC-'Tpr^, (23) 

and 

/'+^ = -KrV//r-, (24) 

respectively. The solution of Eq. (24) is the sum of the solutions of the homo- 
geneous equation 

/' + { = 0, (25) 
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and a particular solution of Eq. (24). These are given by 



/I = — jHT^ (26) 



and 



f2 = -^r^Peffr'. (27) 

The solution /i represents a particle at the origin of coordinates and as such 
is not relevant to our problem. We take, accordingly, /2 as the general solution, 

Using the values of k = SttG/c^t^ and pc = S/SttGt"^, we obtain 

/w = ^-^ (29) 

where flm = p/pc- 

The two solutions given by Eqs. (23) and (29) for /(r) can now be equated, 
giving 

= = = 4.544 (1 - Qm) X 10-'g/cm\ (30) 

Furthermore, from Eqs. (19) and (22) we find that 

e-'^ = 1 + / (r) = 1 + TC-'npr^ = 1 + ^^^r^ (31) 

It will be recalled that the Universe expansion is determined by Eq. (8), 
dr/dv = re"''/^. The only thing that is left to be determined is the signs of 
(1 — flm) or the pressure p. 

Thus we have 



^ = rVl + KTC-'pr^ =r\U + Kfi^r\ (32) 

av V c^T^ 

For simplicity we confine ourselves to the linear approximation, thus Eq. (32) 
yields 

(33) 



dr I ^ n, _i 2 



dv 



T (^1 + — rc pr 



1 I - - -™ 2 

1 H — „ „ „ r 



2c2t2 



5 Classification of universes 

The second term in the square bracket in the above equation represents the 
deviation due to gravity from the standard Hubble law. For without that term, 
Eq. (33) reduces to dr/dv = r, thus r = tv + const. The constant can be taken 
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zero if one assumes, as usual, that at r = the velocity should also vanish. Thus 
r = TV, 01 V = Hqt (since Hq Accordingly, the equation of motion (33) 

describes the expansion of the Universe when 0.^ = 1, namely when p = p^. 
The equation then coincides with the standard Hubble law. 

The equation of motion (33) can easily be integrated exactly by the substi- 
tions 



smx = 
sinh X 



{Qm-l) r 



2cr' 



2 2cT 

One then obtains, using Eqs. (33) and (34), 

dv = cdxl (0„ - 1)^/^ cosx; 

dv = cdxl (1 ^^m)^^^ coshx; 



f^m < 1- 



(34a) 
(346) 

(35a) 
(356) 



We give below the exact solutions for the expansion of the Universe for each 
of the cases, Q,m > 1 and Qm < 1. As will be seen, the case of Clm = 1 can be 
obtained at the limit — ^ 1 from both cases. 

The case flm > !• From Eq. (35a) we have 



/ 



dv 



1)72/ 



dx 

/2J cosx' 



(36) 



where sinx = r/a, and a = ct\J (fi™ — 1) /2. A simple calculation gives [7] 



/ =ln 


1 + sin X 


/ cosx 


cosx 



A straightforward calculation then gives 



V = — In 
2r 



1 + r/a 



1 — r/a 



(37) 



(38) 



As is seen, when r ^ then u ^ and using the L'Hospital lemma, v ^ r/r 
as a ^ (and thus Vim — > 1)- 

The case < !• From Eq. (35b) we now have 



/ 



dv 



dx 



v/(l-0„)/2y coshx' 



(39) 



where sinhx = r/h, and 6 = cr-^/ (1 — fim) /2. A straightforward calculation 
then gives [7] 

.„._x 



/ 



We then obtain 



coshx 
coshx = 



arctan e'* 



1 + 



62' 



(40) 
(41) 
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= sinhx + coshx ~ 5 + V ^ + (^^) 



Equations (39) and (40) now give 

2c 



V(l-J^m)/2 



arctane^ (43) 



where K is an integration constant which is determined by the requirement that 
at r = then v should be zero. We obtain 



K = -7rc/2 V(l - f]„) /2, (44) 

and thus ^ 

V = — ( arctan ) . (45) 

v/(i - n„,) /2 V 4; 

A straightforward calculation then gives 



T 



«=^^2arctan(^ + -^l + ^ (46) 



As for the case Sim > 1 one finds that v ^ when r ^ 0, and again, using 
L'Hospital lemma, r = tv when 6^0 (and thus — * !)• 

6 Physical meaning 

To see the physical meaning of these solutions, however, one does not need the 
exact solutions. Rather, it is enough to write down the solutions in the lowest 
approximation in r~^. One obtains, by differentiating Eq. (33) with respect to 
V, for flm > 1, 

cfr/dv^ = -kr; k = ^"'^J (47) 

the solution of which is 

V V 

r(w)=^sina — hScosa-, (48) 
c c 

where — (fim — 1)/2 and A and B are constants. The latter can be determined 
by the initial condition r (0) = = S and dr {0) /dv = t = Aa/c, thus 

CT V 

r(v) = — sin a-. (49) 
a c 

This is obviously a closed Universe, and presents a decelerating expansion. 
For Qm < 1 we have 

d'r/dv' = ^^J^, (50) 
whose solution, using the same initial conditions, is 

r{v) = jsmh(3-^, (51) 

where = (1 — flm)/^- This is now an open accelerating Universe. 
For Clm = 1 we have, of course, r = tv. 
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7 The accelerating universe 



We finally determine which of the three cases of expansion is the one at present 
epoch of time. To this end wc have to write the solutions (49) and (51) in 
ordinary Hubble's law form v = Hgr. Expanding Eqs. (49) and (51) into power 
series in v/c and keeping terms up to the second order, we obtain 

r = Tv{l- a^v^/Qc^) , (52a) 

r^Tv{l + [3^v^lQc^), {52b) 

for > 1 and < respectively. Using now the expressions for a and /3, 
Eqs. (52) then reduce into the single equation 

r = TV [1 + (1 - Qm) v^l%(?\ . (53) 

Inverting now this equation by writing it as f = Hor, we obtain in the lowest 
approximation 

Ho = h[l-{l- ftm) v^ec^] , (54) 
where h = t~^. To the same approximation one also obtains 

Ho = h[l-{l- flm) z^/&] =h[l-{l- flm) r^/6c^T^] , (55) 

where z is the redshift parameter. As is seen, and it is confirmed by experiments, 
Hq depends on the distance it is being measured; it has physical meaning only at 
the zero-distance limit, namely when measured locally, in which case it becomes 
h = l/T. 

It follows that the measured value of Hq depends on the "short" and "long" 
distance scales [8]. The farther the distance Hq is being measured, the lower 
the value for Hq is obtained. By Eq. (55) this is possible only when Qrn < 1) 
namely when the Universe is accelerating. By Eq. (30) we also find that the 
pressure is positive. 

The possibility that the Universe expansion is accelerating was first predicted 
using CGR by the author in 1996 [2] before the supernovae experiments results 
became known. 

It will be noted that the constant expansion is just a transition stage between 
the decelerating and the accelerating expansions as the Universe evolves toward 
its present situation. 

Figure 1 describes the Hubble diagram of the above solutions for the three 
types of expansion for values of from 100 to 0.245. The figure describes 
the three-phase evolution of the Universe. Curves (l)-(5) represent the stages 
of decelerating expansion according to Eq. (49). As the density of matter p 
decreases, the Universe goes over from the lower curves to the upper ones, 
but it does not have enough time to close up to a big crunch. The Universe 
subsequently goes over to curve (6) with 0,^ = 1, at which time it has a constant 
expansion for a fraction of a second. This then followed by going to the upper 
curves (7) and (8) with Q^n < Ij where the Universe expands with acceleration 



10 



according to Eq. (51). Curve no. 8 fits the present situation of the Universe. 
For curves (l)-(4) in the diagram we use the cutoff when the curves were at 
their maximum. In Table 1 we present the cosmic times with respect to the big 
bang, the cosmic radiation temperature and the pressure for each of the curves 
in Fig. 1. 

Figures 2 and 3 show the Hubble diagrams for the distance-redshift rela- 
tionship predicted by the theory for the accelerating expanding Universe at the 
present time, and Figures 4 and 5 show the experimental results. 

Our estimate for h, based on published data, is /i w 80 km/sec- Mpc. As- 
suming « 80 km/sec-Mpc, Eq. (55) then gives 

Ho = h[l- 1.3 X lO"'' (1 - 0„) r^] , (56) 

where r is in Mpc. A computer best-fit can then fix both h and ^m- 

To summarize, a theory of cosmology has been presented in which the dy- 
namical variables are those of Hubble, i.e. distances and velocities. The theory 
descirbcs the Universe as having a three-phase evolution with a decelerating ex- 
pansion, followed by a constant and an accelerating expansion, and it predicts 
that the Universe is now in the latter phase. As the density of matter decreases, 
while the Universe is at the decelerating phase, it does not have cinough time to 
close up to a big crunch. Rather, it goes to the constant-expansion phase, and 
then to the accelerating stage. As we have seen, the equation obtained for the 
Universe expansion, Eq. (51), is very simple. 



8 Theory versus experiment 

The Einstein gravitational field equations with the added cosmological term are 
[9]: ^ 

R^iu - -j^giivR + A.g,„y = kT^^, (57) 

where A is the cosmological constant, the value of which is supposed to be 
determined by experiment. In Eq. (57) i?^^ and R are the Ricci tensor and 
scalar, respectively, k = 87rG, where G is Newton's constant and the speed of 
light is taken as unity. 

Recently the two groiips (the Supernovae Cosmology Project and the High-Z 
Supernova Team) concluded that the expansion of the Universe is accelerating 
[10-16]. The two groups had discovered and measured moderately high redshift 
(0.3 < z < 0.9) supernovae, and found that they were fainter than what one 
would expect them to be if the cosmos expansion were slowing down or constant. 
Both teams obtained 

« 0.3, r!A » 0.7, (58) 

and ruled out the traditional {ilm, ^^a) = (1, 0) Universe. Their value of the 
density parameter Oa corresponds to a cosmological constant that is small but, 
nevertheless, nonzero and positive, 

A w lO-^^m-^ w IQ-^h-^. (59) 
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In previous sections a four-dimensional cosmological theory (CGR) was pre- 
sented. Although the theory has no cosmological constant, it predicts that the 
Universe accelerates and hence it has the equivalence of a positive cosmological 
constant in Einstein's general relativity. In the framework of this theory (see 
Section 2) the zero-zero component of the field equations (3) is written as 

R°0-ls°oR = ^Peff = K{p-pc), (60) 

where pc = ^/kt^ is the critical mass density and r is Hubble's time in the 
zero-gravity limit. 

Comparing Eq. (60) with the zero-zero component of Eq. (57), one obtains 
the expression for the cosmological constant of general relativity, 

A = Kpc = (61) 

To find out the numerical value of t wc use the relationship between h = 
and Hq given by Eq. (55) (CR denote values according to Cosmological 
Relativity): 

Ho = h[l-{l-ilg^)zy6], (62) 

where z = v/c is the rcdshift and fJ^^ = Pm/Pc with pc = Sh'^/SirG. (Notice 
that our pc = 1.194 x 10~^^ g/cm^ is different from the standard pc defined 
with Hq.) The redshift parameter z determines the distance at which Hq is 
measured. We choose z = 1 and take for 

ng"" = 0.245, (63) 

its value at the present time (see Table 1) (corresponds to 0.32 in the standard 
theory), Eq. (62) then gives 

Ha = 0.874ft,. (64) 

At the value z = 1 the corresponding Hubble parameter according to the 
latest results from HST can be taken [17] as Hq = 70km/s-Mpc, thus h = 
(70/0.874)km/s-Mpc, or 

h = 80.092km/s-Mpc, (65) 

and 

T = 12.486Gyr = 3.938 x lO^^s. (66) 

What is left is to find the value of ^1^- We have — pf^ / p^, where 
p^J = SH^/SttG and p^ = Sh'^/SirG. Thus = {Ho/hf = 0.874^, or 

OC'ii = 0.764. (67) 
As is seen from Eqs. (63) and (67) one has 

nr = + r^A^ = 0.245 + 0.764 = 1.009 w 1, (68) 
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which means the Universe is Euchdean. 

As a final result we calculate the cosmological constant according to Eq. 
(61). One obtains 

A = S/r^ = 1.934 x IQ-^^s'^. (69) 

Our results confirm those of the supernovae experiments and indicate on 
the existance of the dark energy as has recently received confirmation from the 
Boomerang cosmic microwave background experiment [18,19], which showed 
that the Universe is Euclidean. 

9 Some remarks 

In this paper the cosmological general relativity, a relativistic theory in spaceve- 
locity, has been presented and applied to the problem of the expansion of the 
Universe. The theory, which predicts a positive pressure for the Universe now, 
describes the Universe as having a three-phase evolution: decelerating, constant 
and accelerating expansion, but it is now in the latter stage. Furthermore, the 
cosmological constant that was extracted from the theory agrees with the ex- 
perimental result. Finally, it has also been shown that the three-dimensional 
spatial space of the Universe is Euclidean, again in agreement with observations. 

Recently [20,21], more confirmation to the Universe accelerating expansion 
came from the most distant supernova, SN 1997ff, that was recorded by the 
Hubble Space Telescope. As has been pointed out before, if we look back far 
enough, we should find a decelerating expansion (curves 1-5 in Figure 1). Beyond 
z = I one should see an earlier time when the mass density was dominant. The 
measurements obtained from SN 1997ff's redshift and brightness provide a direct 
proof for the transition from past decelerating to present accelerating expansion 
(see Figures 6 and 7). The measurements also exclude the possibility that the 
acceleration of the Universe is not real but is due to other astrophysical effects 
such as dust. 

Table 2 gives some of the cosmological parameters obtained here and in the 
standard theory. 

10 Comparison with general relativity 

In order to compare the present theory with general relativity, we now add 
the time coordinate. We then have a time-space-velocity Universe with two 

time-like and three space-like coordinates, with signature (H h). We 

will be concerned with the classical experiments of general relativity and the 
gravitational waves predicted by that theory. In the following we show that all 
these results are also obtained from the present theory. To this end we proceed 
as follows. 

We first find the cosmological-equivalent of the Schwarzschild spherically- 
symmetric solution in cosmology. It will be useful to change variables from the 
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classical Schwarzschild metric to new variables as follows: 



sin^ X = '>'s/i', dr = —Ivg sin ^ xcosx^X, 



(70) 



where Vg = 2GM/(? is the Schwarzschild radius. We also change the time 
coordinate cdt = Vsdr], thus 77 is a time parameter. The classical Schwarzschild 
solution will thus have the following form: 

ds'^ = [cos^ xdr/2 - 4 sin"*^ xdx^ - sin"^ x [dO"^ + sin^ Odcj)^)] . (71) 

So far this is just the classical spherically symmetric solution of the Einstein field 
equations in four dimensions, though written in new variables. The non-zero 
Christoffel symbols are given by 

Tqi = -sinxcos"^x, ^00 = - ;| sin'' X cos x, 



-3 sin xcosx, 



^ 22 



-sm xcosx, 



^ 33 



- sin X cos X sin 6, 



p2 

^ 12 



-2 sin xcosx, 



(72) 



p2 

J- 33 



■ Sin y cos t 



p3 
J- 13 



-2 sin xcosx. 



r^3 = sin ^9 cos 9. 



It is very lengthy, but one can verify that all components of the Ricci tensor 
i?a/3 are equal to zero identically. 

We now extend this solution to cosmology. In order to conform with the 
standard notation, the zero component will be chosen as the time parameter, 
followed by the three space-like coordinates and then the fourth coordinate 
representing the velocity rdv. We will make one more change by choosing 
rdv = rgdu, thus u is the velocity parameter. The simplest way to have a 
cosmological solution of the Einstein field equation is using the so-called co- 
moving coordinates in which: 

rfs^ = [cos^ xdv^ - 4 sin"*^ xdx^ - sin"'' X {dO^ + sin^ + du^] . (73) 

The coordinates are now x° = 7], = x, =6, = (j), and = u, and Tg is 

now a function of the velocity u, = (u) to be determined by the Einstein 

field equations in five dimensions. Accordingly we have the following form for 
the metric: 



I cos^x 



V 



-4sin-^- 



sin 



— sin X sin 9 



\ 



1 / 



(74a) 



-g = 2rf sin x cos x sin 9. 



(746) 
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The non-zero Christoffel symbols are given by 
Tqi = - sinxcos"^ x, ^04 = 



00 



--sm xcosx, 



^ 11 



-3 sin xcosx, 



^ 14 



^ 22 



1 

-sm XCOSX, 



^ 33 



1 . .2 

-Sin XCOSX Sin 



p2 
^ 12 



2 sin xcosx, 
rfa = -2sin"^xcosx, 



p2 _ ■ -1 
1 24 — ' s's , 



33 



sm tf cos f 



(75) 



^ 00 



-r^r^ ^ cos^ X, 



11 



r^3 = sin ^6 cos 9, 
-- 4r«r7^ sin"^x, 



p3 
^ 34 



22 



• -1 • -4 

r=r„ sm 



X, 



r33 = rsr, ^ si 



sin X sin^ ( 



where the dots denote derivatives with respect to the velocity parameter u. 

The Ricci tensor components after a lengthy but straightforward calculation, 
are given by: 



-Roo = - {rsTs ' + 2rs^r^ ^) cos^ x, 



Rii = 4 {rSr-^ + 2r; 
-R22 = {r'sr, 
R33 = 



^r^^)sm ^x, 



2r;V^2)sin ^ x. 



(76) 



(r'sTj, ^ + 2r's^rg ^) sin ^ x sin^ 



-4(r 



• 2 —2 

Ts r. ^ 



)• 



All other components are identically zero. 

We are interested in vacuum solution of the Einstein field equations for 
the spherically symmetric metric (generalized Schwarzschild to cosmology) , the 
right-hand sides of the above equations should be taken zero. A simple calcula- 
tion then shows that = 0, =0. Accordingly the cosmological Schwarzschild 
metric is given by Eq. (74a) with a constant Vg = 2GM/c^. The metric (74a) 
can then be written, using the coordinate transformations (70), as 



I 

r 



V 



(1-1^)" 



— sin^ I 



\ 



1 / 



(77) 



r, x"^ = 9, = 



= 6, and 



where the coordinates are now xP = ct, x^ = 

We are now in a position to compare the present theory with general rela- 
tivity. 
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11 Gravitational redshift 



We start with the simplest experiment of the gravitational redshift. Although 
this experiment is not considered as one of the proofs of general relativity (it 
can be derived from conservation laws and Newtonian theory). 

Consider two clocks at rest at two points denoted by 1 and 2. The prop- 
agation of light is detc^rniined by ds at each point. Since at these points all 
spatial infinitesimal displacements and change in velocities vanish, one has 
ds^ = gooc^dt^. Hence at the two points we have 

dsil)^[goo{l)f^cdt, (78a) 

ds (2) = [goo {2)f^cdt (786) 

for the proper time (see Fig. 8). 

The ratio of the rates of similar clocks, located at different places in a grav- 
itational field, is therefore given by 

ds (2) /ds (1) = [.goo (2) /goo (1)]'/' . (79) 

The frequency uo of an atom located at point 1, when measured by an observer 
located at point 2, is therefore given by 

^^ = ^^o[5oo(l)/5oo(2)]'/'. (80) 

If the gravitational field is produced by a spherically symmetric mass distri- 
bution, then we may use the generalized Schwarzschild metric given above to 
calculate the above ratio at the two points. In this case 500 = 1 — 2GM/(?r, 
and therefore 

[500 (1) /500 (2)]'/' « 1 + (GM/c^) (l/r2 - 1/ri) 
to first order in GM/(?r. We thus obtain 

Aiz/z/o ^{v- vo) /vo ~ - {GM/c") (l/n - 1/^2) 

for the frequency shift per unit frequency. Taking now ri to be the observed 
radius of the Sun and r2 the radius of the Earth's orbit around the Sun, then 
we find that 

Aiy/i^o « -GMsun/c^rsun, (81) 

where Msun and rsun are the mass and radius of the Sun. Accordingly we 
obtain Av/vq « —2.12 x 10~^ for the frequency shift per unit frequency of the 
light emitted from the Sun. The calculation made above amounts to neglecting 
completely the Earth's gravitational field. The above result is the standard 
gravitational redshift (also known as the gravitational time dilation). 
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12 Motion in a centrally symmetric gravitational 
field 



We assume that small test particles move along geodesies in the gravitational 
field. We also assume that planets have small masses as compared with the mass 
of the Sun, to the extent that they can be considered as test particles moving in 
the gravitational field of the Sun. As a result of these assumptions, the geodesic 
equation in the cosmological Schwarzschild field will be taken to describe the 
equation of motion of a planet moving in the gravitational field of the Sun. In 
fact, we do not need the exact solution of the generalized Schwarzschild metric 
(77), but just its first approximation. We obtain in the first approximation the 
following expressions for the components of the metric tensor; 

goo = l-rs/r, gom = 0, 9o4=0, 

Qmn = -Smn " rsX™x"/r^, .g„4 = 0, 044 = 1- (82a) 

The contravariant components of the metric tensor are consequently given, in 
the same approximation, by 

g°° = l + r,/r, 5°" = 0, 5°^ = 0, 

,g™" = -5™" + r,x"x"/r3, c,™^ ^ ^44 ^ ^ ^^^b) 

We may indeed verify that the relation g^ixg^" = between the contravariant 
and covariant components of the above approximate metric tensor is satisfied to 
orders of magnitude of the square of rg/r. A straightforward calculation then 
gives the following expressions for the ChristofTel symbols: 

_ d (\ 

■i- — 



°" 2 9a;" V r 



^mn = 1's^Omn - -jTg ^ . 

All other components vanish. 

We now use these expressions for the Christoffel symbols in the geodesic 
equation 

+ - T%xk) x-xf^ = 0, (84) 

where a dot denotes differentiation with respect to the time coordinate x^. We 
obtain 
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110 

r, d (\ 



2 dx'' 



ix-'x'')^ x^ 



( - 1 - i-^x^ — (-)- — 

2r dx^ \r J dx'^ \r J 2r^ 

(856) 

Consequently we obtain from the geodesic equation (84) the following equation 
of motion for the planet: 

••^ rs d f V 



^ Sx*^ i r / 2r dx^ \ r) ^ \ r) ^ ^ 2r'° ("^ ^ ) 



a; 



(86) 



Replacing now the derivatives with respect to x° by those with respect to f(E 
x°/c) in the latter equation, we obtain 



X - GMV- = r, 

r 



(87) 

where use has been made of the three-dimensional notation. 

Hence the equation of motion of the planet differs from the Newtonian one 
since the left-hand side of Eq. (87) is proportional to terms of order of magnitude 

instead of vanishing identically. This correction leads to a fundamental effect, 
namely, to a systematically secular change in the perihelion of the orbit of the 
planet. 

To integrate the equation of motion (87) we multiply it vectorially by the 
radius vector x. We obtain 

xxx = -rs(x- V(l/r))(xxx). (88) 

All other terms in Eq. (87) are proportional to the radius vector x and thus 
contribute nothing. Equation (88) may be integrated to yield the first integral 

xxx = Je-'■'/^ (89) 

Here J is a constant vector, the angular momentum per mass unit of the planet. 
One can easily check that the first integral (89) indeed leads back to Eq. (88) 
by taking the time derivatives of both sides of Eq. (89). 

Prom Eq. (89) we see that the radius vector x moves in a plane perpendicular 
to the constant angular momentum vector J, thus the planet moves in a plane 
similar to the case in Newtonian mechanics. If we now introduce in this plane 
coordinates r and cj) to describe the motion of the planet, the equation of motion 
(87) consequently decomposes into two equations. Introducing now the new 
variable u = 1/r, we can then rewrite the equations in terms of m(^), using 

r = KfP, 



12 



2u^ -2 ^* 12 " 
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where a prime denotes a differentiation with respect to the angle (j). We subse- 
quently obtain 

•• u' -2 2GM , -2 

u & 

A straightforward calculation then gives, using the expression for 

,2^ 2' 



2 

u - GM ( ^ ' 



u^Y GM 



2^2 _ y'2 _ 2GMu [ ^ 



(90). 



The latter equation can be further simplified if we use the first integral 

r2^ = Je-2GM/-^-. 



We obtain 



6 J 



■ ' ~ ~ J2 " 



Hence, to an accuracy of 1/c , Eq. (90) gives 

GM GM ( 2 ,2 GM \ , , 

« + " - ^ = ^ (2«^ - « + 2^"j • (91) 

Equation (91) can be used to determine the motion of the planet. The 
Newtonian equation of motion that corresponds to Eq. (91) is one whose left- 
hand side is identical to the above equation, but is equal to zero rather than to 
the terms on the right-hand side. This fact can easily be seen if one lets GM/c? 
go to zero in Eq. (91). Therefore in the Newtonian limit we have 

«" + «-^«0, (92) 
whose solution can be written as 

uq(\ + e cos (j)) . (93) 

Here Mo is a constant, and e is the eccentricity of the ellipse, e = (1 — b"^ /a?)^/"^, 
where a and h are the semimajor and semiminor axes of the ellipse. Using 
the solution (93) in the Newtonian limit of the equation of motion (92) then 
determines the value of the constant mq, as uq = GM/ ,P . 

To solve the equation of motion (91), we therefore assume a solution of the 
form 

u = uo{l + e cos acj)) , (94) 

where a is some parameter to be determined, and whose value in the usual 
nonrelativistic mechanics is unity. The appearance of the parameter a 1 in 
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our solution is an indication that the motion of the planet will no longer be a 
closed ellipse. 

Using the above solution in Eq. (91), and equating coefficients of cos acf), 
then gives 



0^ = 1 P- ( 2uo + 



2GM ( GM 



j2 ; ■ 

If we substitute for GM/ J'^ in the above equation its nonrelativistic value uq, 
then the error will be of a higher order. Hence the latter equation can be written 
as 

a =1 ^tto 

or 

, 3GM 

a = 1 5— Wo- (95) 

Successive perihelia occur at two angles (/>! and ^2 when a<f)2 — ctipi = 27r. 
Since the parameter a is smaller than unity, we have 4>2 — (j)! — 27r/a > 2tt. 
Hence we can write (/)2 — 0i = 27r + A(j), with Acj) > 0, or 

(3GM \ 
1 —Uo J (27r + Acj)) = 2w. (96) 

As a result there will be an advance in the perihelion of the orbit of the planet 
per revolution given by Eq. (96) or, to first order, by 

A(/> = 67rGMuo/c^ (97) 

The constant uq can also be expressed in terms of the eccentricity, using the 
Newtonian approximation. Denoting the radial distances of the orbit, which 
correspond to the angles <p2 = and (pi = w, by r2 and ri , respectively, we have 
fromEq. (93), 

l/r2 = Uo (1 + e) , 1/ri = Uo (1 - e) . 
Hence since ri + r2 = 2a, we obtain (see Fig. 9) 

2a = n + r2 = 2/wo (l - e^) , 

where a is the semimajor axis of the orbit, and therefore 

Uo = 1/a (1 — e^) . 

Using this value for mq in the expression (97) for A(f), we obtain for the perihelion 
advance the expression 

c^a (1 — e^) 

in radians per revolution (see Fig. 10). This is the standard general relativistic 
formula for the advance of the perihelion. 

In the next section we discuss the deflection of a light ray moving in a 
gravitational field. 



20 



13 Deflection of light in a gravitational fleld 

To discuss the effect of gravitation on the propagation of hght signals we may use 
the geodesic equation, along with the null condition ds = at a fixed velocity. 
A light signal propagating in the gravitational field of the Sun, for instance, will 
thus be described by the null geodesies in the cosmological Schwarzschild field 

at dv = 0. 

Using the approximate solution for the cosmological Schwarzschild metric, 
given by Eq. (82a), we obtain 



Hence we have, to the first approximation in GM/c^, the following equation of 
motion for the propagation of light in a gravitational field: 



where a dot denotes differentiation with respect to the time coordinate t{= 



Just as in the case of planetary motion (see previous section), the motion 
here also takes place in a plane. Hence in this plane we may introduce the polar 
coordinates r and ^. The equation of motion (100) then yields, to the first 
approximation in GM/c^, the following equation in the polar coordinates: 




= 0. (99) 




(100) 




(101) 



Changing now variables from r to u{(j)) = 1/r, we obtain 




(102) 



where a prime denotes differentiation with respectto the angle (p. 
Moreover we may use the first integral of the motion, 



(103) 



in Eq. (102), thus getting 




(104) 



Differentiation of this equation with respect to </> then gives 



u" + U + 



2GM 



(105) 
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In Eq. (105) terms have been kept to the first approximation in GMjc^ only. 

To solve Eq. (105) we notice that, in the lowest approximation, we have, 
fromEq. (104), 

u" w -u. (107) 
Hence using these approximate expressions in Eq. (105) gives 

m" + « = — 5— (108) 

for the equation of motion of the orbit of the light ray propagating in a spheri- 
cally symmetric gravitational field. 

In the lowest approximation, namely, when the gravitational field of the 
central body is completely neglected, the right-hand side of Eq. (108) can be 
taken as zero, and therefore u satisfies the equation u" -\-u = 0. The solution 
of this equation is a straight line given by 

u=\sm4>, (109) 

where i? is a constant. This equation for the straight lino shows that r = l/u 
has a minimum value R at the angle ^ = 7r/2. If we denote y = rsin^, the 
straight line (109) can then be described by 

y = r sin (/) = ii = constant (HO) 

(see Fig. 11). 

We now use the approximate value for u, Eq. (109), in the right-hand side of 
Eq. (108), since the error introduced in doing so is of higher order. Wo therefore 
obtain the following for the equation of motion of the orbit of the light ray: 

„ 3GM . 2 , 

(111) 

The solution of this equation is then given by 

u= -^sin<^+^ (l + cos2<^) . (112) 

Introducing now the Cartesian coordinates x = rcos^ and y = rsin0, the 
above solution can then be written as 

GM 2x2+^2 

We thus see that for large values of |x| the above solution asymptotically ap- 
proaches the following expression: 

y«i?-^N. (114) 
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As seen from Eq. (114), asymptotically, the orbit of the light ray is described 
by two straight lines in the spacetime. These straight lines make angles with 
respect to the x axis given by tan^ = ± (2GM/c^i?) (see Fig. 12). The angle 
of deflection A(f) between the two asymptotes is therefore given by 

A^=^. (115) 

This is the angle of deflection of a light ray in passing through the grav- 
itational field of a central body, described by the cosmological Schwarzschild 
metric. For a light ray just grazing the Sun, Eq. (115) gives the value 

AS = = 1.75seconds. 

c'-^Sun 

This is the standard general-relativistic formula. Observations indeed confirm 
this result. One of the latest measurements gives 1.75±0.10 seconds. It is worth 
mentioning that only general relativity theory and the present theory predict 
the correct factor of the deflection of light in the gravitational field. 

In the next section the gravitational radiation prediction is considered. 



14 Gravitational radiation 

In the following we show that the present theory also predicts gravitational 
radiation, a distinguished result of classical general relativity theory. We will 
not develop a complete theory of gravitational radiation. Rather we will confine 
ourselves in showing that the present theory does predict the phenomenon. This 
is done in the weak field approximation, as is usually done in standard general 
relativity theory. 

14.1 Linear approximation 

For convenience, the coordinate system to be used in the linearized theory will 
be Cartesian, and hence the Minkowskian metric will have the form 

7?^,= 77^-^ = (1,-1, -1,-1,1), (116) 

when c and t are taken as unity. The gravitational field described by the metric 
tensor g^i, is now called weak if it differs from the Minkowskian metric tensor 
by terms which are much smaller than unity, 

\9^^u - ri^A < 1. (117) 

The above condition need not be satisfied in the entire spacetime, and it could 
be valid at a region of it. 

We now assume that the metric tensor can be expanded as an infinite series, 

9nv = Vniy + >^\9nv + >^^29nv H , (118) 
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where A is some small parameter, and we limit ourselves to the first-order term 
alone. Hence we can write 

~ ViJ-i' + h^v, (119a) 

where h^v = ^iQnv We also expand the contravariant components of the metric 
tensor, 

g^"" K -q^"" + h'"' . (11%) 

Prom the condition g^xg^" = S'^ one then is able to relate /i'*'' to h^i, (neglecting 
nonlinear terms), 

hf"" = -r]'"'r]'"'Kp. (120) 
14.2 The linecirized Einstein equations 

We can now derive the linearized Einstein equations. To this end we have to 

find the first approximate value of the Einstein tensor, the Ricci tensor, the 
Ricci scalar, and the Christoffel symbols. A simple calculation then gives 

Kfi ~ Iv"^ ihxa,0 + hxp,a - K0,X) (121) 
for the Christoffel symbols and 

Raises ~ 2 i^o:S,0-f + hffj^aS — ^/35,a7 — ha^^ffs) (122) 

for the Riemann tensor. Accordingly we have the following expressions for the 
Ricci tensor, the Ricci scalar, and the Einstein tensor, respectively: 



^IJ.iy ^ 2 [^'/^^ ^^l^^^P<7 hfxp-jUa ^i/p,/x(t) Vfii'V^ (^^^pcr V ^ ^pcr,afs')'\ ; 



R « {hc,5,0'y - hi3S,a^) (124) 

(125) 
where h = rj^^haj^- 

A simplification in the linearized field equations occurs if we introduce the 
new variables ^ 

7p,i^ = hfj,^ - -rjij.t^h, (126) 

from which one obtains ^ 

h^,^ = -Jpu - ^Vp-ul, (127) 

with 7 = T]"'^jap- Introducing the above expressions into the Einstein field 
equations we obtain 

07m. - ^"'^ (7aM,/3^ + 7a^,/3^) + Vf..^V^''v"^7Xa,p0 = -2kT^. (128) 
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for the linearized gravitational field equations. In Eq. (128) the symbol Q is 
the operator in flat space, 

0/ = rt..= (^^-V^ + ^|^)/. (129) 

We can simplify still further the above field equations by choosing coordi- 
nates in which 

7^ = V'"''r^.p,. = 0. (130) 

This is similar to choosing a gauge in solving the wave equation in electrody- 
namics. As a result we finally obtain for the linearized Einstein equations the 
following: 

07m. = (131) 
along with the supplementary condition 

V'^I^.P,. = 0, (132) 

which solutions 7^1^ of Eq. (131) should satisfy. Finally we see from Eq. (131) 
that a necessary condition for Eq. (132) to be satisfied is that 

v'^'^Tf,^,!} = 0, (133) 

which is an expression for the conservation of the energy and momentum without 
including gravitation. 



14.3 Gravitational waves 
In vacuum, Eq. (131) reduces to 

07m. = 0, (134) 



or 



1 52 



c 



Thus the gravitational field, like the electromagnetic field, propagates in vac- 
uum with the speed of light. The above analysis also shows the existance of 
gravitational waves. 
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Tabic 1: The Cosmic Times with respect to the Big Bang, the Cosmic Tem- 
perature and the Cosmic Pressure for each of the Curves in Fig. 1. 



Curve Clm Time in Units Time Temperature Pressure 
No*. ofr (Gyr) (K) {g/cm^) 







DECELERATING EXPANSION 




1 


100 


3.1 X IQ-*^ 


3.87 X 10-5 


1096 


-4.499 


2 


25 


9.8 X 10-5 


1.22 X 10-3 


195.0 


-1.091 


3 


10 


3.0 X 10-4 


3.75 X 10-3 


111.5 


-0.409 


4 


5 


1.2 X 10-3 


1.50 X 10-2 


58.20 


-0.182 


5 


1.5 


1.3 X 10-2 


1.62 X 10-1 


16.43 


-0.023 






CONSTANT EXPANSION 






6 


1 


3.0 X 10-2 


3.75 X 10-1 


11.15 









ACCELERATING EXPANSION 




7 


0.5 


1.3 X 10-1 


1.62 


5.538 


+0.023 


8 


0.245 


1.0 


12.50 


2.730 


-h0.034 



*The calculations are made using Carmcli's cosmological transformation, 
Eq. (2), that relates physical quantities at different cosmic times when gravity 
is extremely weak. 

For example, we denote the temperature by 6, and the temperature at the 
present time by 9o, we then have 

Q _ Oq _ 00 _ 2.7SK _ 2.73K 

where T is the time with respect to B.B. 

The formula for the pressure is given by Eq. (30), p = c{l — Q) / SttGt . Using 
c = 3 X 10i°cm/s, r = 3.938 x lO^'^s and G = 6.67 x IQ-^cm^/gs^, we obtain 

p = 4.544 X 10-2 (1 - Ci) g/cm^. 



28 



Table 2: Cosmological parameters in cosmological general relativity and in 

standard theory 





COSMOLOGICAL 


STANDARD 




RELATIVITY 


THEORY 


Theory type 


Spacevelocity 


Spacetime 


Expansion 


Tri-phase: 


One phase 


type 


decelerating, constant, 
accelerating 




Present expansion 


Accelerating 


One of three 




(predicted) 
0.034g/cm2 


possibilities 


Pressure 


Negative 


Cosmological constant 


1.934 X 10-35s-2 
(predicted) 


Depends 


r^T^ = ^jYi -\- 


1.009 


Depends 


Constant-expansion 


8.5Gyr ago 


No prediction 


occurs at 






Constant-expansion 


Fraction of 


Not known 


duration 


second 




Temperature at 


146K 


No prediction 


constant expansion 
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FIGURE CAPTIONS 



Fig. 1 Hubblc's diagram describing the three-phase evolution of the Universe 
according to cosmological general relativity theory. Curves (1) to (5) represent 
the stages of decelerating expansion according to r{v) = (cr/a) sin av/c, where 
= (n — l)/2, fl = p/pc with pc a constant, pc = S/SttGt^, and c and r are 
the speed of light and the Hubble time in vacuum (both universal constants). 
As the density of matter p decreases, the Universe goes over from the lower 
curves to the upper ones, but it does not have enough time to close up to a 
big crunch. The Universe subsequently goes to curve (6) with = 1, at which 
time it has a constant expansion for a fraction of a second. This then followed 
by going to the upper curves (7), (8) with f7 < 1 where the Universe expands 
with acceleration according to r{v) = {ct/(3) smh(}v/c, where = (1 — f2)/2. 
Curve no. 8 fits the present situation of the Universe. (Source: S. Behar and 
M. Carmeli, Ref. 3) 

Fig. 2 Bubble's diagram of the Universe at the present phase of evolution with 
accelerating expansion. (Source: S. Behar and M. Carmeli, Ref. 3) 

Fig. 3 Hubblc's diagram describing decelerating, constant and accelerating ex- 
pansions in a logarithmic scale. (Source: S. Behar and M. Carmeli, Ref. 3) 

Fig. 4 Distance vs. redshift diagram showing the deviation from a constant 
toward an accelerating expansion. (Source: A. Riess et al, Ref. 12) 

Fig. 5 Relative intensity of light and relative distance vs. redshift. (Source: A. 
Riess et al., Ref. 12) 

Fig. 6 Hubble diagram of SNe la minus an empty (i.e., "empty" f2 = 0) Uni- 
verse compared to cosmological and astrophysical models. The points are the 
redshift-binned data from the HZT (Riess et al. 1998) and the SCP (Perlmutter 
et al. 1999). The measurements of SN 1997ff are inconsistent with astrophysical 
effects which could mimic previous evidence for an accelerating Universe from 
SNe la at z « 0.5. (Source: A. Riess et al, Ref. 21) 

Fig. 7 Same as Fig. 6 with the inclusion of a family of plausible, flat f^A 
cosmologies. The transition redshift (i.e., the coasting point) between the ac- 
celerating and decelerating phases is indicated and is given as [2flA/^MY^^ — 1- 
SN 1997ff is seen to lie within the epoch of deceleration. This conclusion is 
drawn from the result that the apparent brightness of SN 1997ff is inconsistent 
with values of JIa > 0.9 and hence a transition redshift greater than that of SN 
1997ff. (Source: A. Riess et al, Ref. 21) 

Fig. 8 Propagation of light in curved spacetime. 

Fig. 9 Newtonian limit of planetary motion. The motion is described by a 
closed ellipse if the effect of other planets is completely neglected. 

Fig. 10 Planetary elliptic orbit with perihelion advance. The effect is a general 
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relativistic one. The advance of the perihelion is given by A(/) in radians per 
revolution, where A(/> = QttGM / c^a{l — e^), with M being the mass of the Sun, 
a the semimajor axis, and e the eccentricity of the orbit of the planet. 

Fig. 11 Light ray when the effect of the central body's gravitational field 
is completely neglected. The light ray then moves along the straight line 
y = rsincj) = R = constant, namely, u = 1/r = (1/i?) sine/; 

Fig. 12 Bending of a light ray in the gravitational field of a spherically sym- 
metric body. The angle of deflection At/) = AGM/c^R, where M is the mass of 
the central body and R is the closest distance of the light ray from the center 
of the body. 
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Planet 
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